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Abstract , hagnetospheric studies often require knowledge of the orien- 
tation of the IMF. In order to test the accuracy of using magnetometer 
data from a spacecraft orbiting the sunward libratlon point for this 
purpose, the angle between the IMF at ISEE 3, when It was positioned around 
the libratlon point, and at ISEE 1, orbiting Earth, has been calculated for 
a data set of two-hour periods covering four months. For each period, a 
ten-minute average of ISEE 1 data is compared with ten-minute averages of 
ISEE 3 data at successively lagged intervals. At the lag time equal to the 
time required for the solar wind to convect from ISEE 3 to ISEE 1, the 
median angle between the IMF orientation at the two spacecraft is 20°, end 
80% of the cases have Bngles less than 38°. The results for the angles 
projected on the y-z plane are essentially the same. The minimum angle 
between the IMF orientation at the two spacecraft has a median of 11° , with 
80% less than 19°. These low values indicate little temporal or small scale 
variation between the spacecraft. The minimum angle generally occurs at a 
lag time different from the convection time. The sign of the difference 
depends on IMF orientation in the sense that magnetic features tend to 
arrive sooner when the IMF is directed along the line between the space- 
craft. However, the difference between a lag time appropriate to this 
corotation geometry and the convection lag time is not large enough to 
produce a significant decrease in the angles between the IMF vectors at the 
two spacecraft. We conclude that the IMF at a libration-point-orblting 
spacecraft, lagged by the time required for the solar wind to convect to 
the earth, is a good, convenient predictor of the IMF near the earth. 


Introduction 


Knowledge of the orientation of the interplanetary magnetic field (IMF) 
Just upstream of the earth is essential for studies of its interaction with 
the magnetosphere [e.g., Burch, 1983]. In the past, satellites in earth 
orbit could monitor the IMF only along the limited portions of their orbits 
which fell outside the bow shock. Continuous measurements of the IMF became 
available with the launch ISEE 3 in 1978 into a halo orbit around the libra- 
tion point, "230 R^. upstream. However, from this position the Bolar wind 
requires about one hour to convect downstream to the earth. Since the IMF 
is’ highly variable In time and space, it is Important to test how accurately 
the ISEE 3 measurements can represent the IMF in earth's vicinity. 

Two studies of the correlations between magnetometer measurements from 
ISEE 1 and ISEE 3 already have been performed. Russell et al. L19803 found 
that the degree of correlation often changes drastically in consecutive 
three-hour intervals, and Crooker et al. 119821 show that the degree of cor- 
relation Increases with the variance of the IMF and decreases significantly 
for spacecraft separations perpendicular to the earth-sun line greater than 
90 Rjr. Also, King 11983] found relatively good correlations between hourly 
averages of measurements made by the earth-orbiting IMP 8 and ISEE 3 for a 
wide range of solar wind parameters. 

The present study addresses more directly the question of how well the 
IMF orientation can be predicted near the earth from the upstream measure- 
ments. The angle between field orientations at ISEE 1 and ISEE 3 is calcu- 
lated for two predicted lag times and compared to the minimum angle between 
the fields in the two-hour interval following the time of the ISEE 3 meas- 
urement. The minimum angle gives the limit on how well the IMF orientation 


can be predicted at any given time. HiBtograms of these angles ore then 
compiled from measurements covering a four-month period, and from them one 
can obtain the probability that the IMF orientation at the upstream libra- 
tion point is within a certain angle of its orientation near the earth. 
Also statistical Information is obtained on the accuracy of predicted time 
lags. 


Data 

Nearly continuous magnetometer data from ISEE 1 and ISEE 3 were 
obtained for the period from July 25 to Nov. 27, 1979. Times chosen for 
analyses were limited to intervals when the earth orbiting ISEE 1 was 
located upstream of the bow shock. The spatial coverage of the satellite 
orbits during this period are shown in Figure 1 in GSE coordinates. The 
orbit of ISEE 3 ranged from 195 to 265 R E along the x axis and reached 
distances of approximately 100 R E away from the x axis, in the y-z plane. 
The x distances correspond to solar wind transit times to earth on the 
order of one hour for typical solar wind speeds. The shaded region indi- 
cates coverage of the ISEE 1 orbit. It shows that the intersatellite 
separation distances in the y-z plane ranged from 18 to 10S Rj?. 

Measurements of solar wind velocity were required to calculate the 
solar wind transit times from ISEE 3 to ISEE 1. Ninety— minute averages of 
velocity were obtained from the Los Alamos solar wind experiment on ISEE 3 
(supplied to data pool by S. J. Bane, principal investigator). 
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Analysis and Results 

Time Variations of Angular Differences . The angle between the magnetic 
field orientations at ISEE 1 and ISEE 3 was calculated as a function of lag 
time for 80S two-hour data segments overlapped by one hour. A ten-minute 
average of ISEE 3 data was compared to successive ten-minute running 
averages of ISEE 1 data at one-minute intervals. Angles were calculated 
both between the total magnetic field vectors and between their projections 
in the y-z plane, perpendicular to the earth-sun line. 

Averages over ten minutes were chosen because ten minutes is roughly 
the minimum length of time required for the magnetosphere to respond to 
changes in IMF orientation. Averages over S, 20, and 60 min were tested; 
but the shorter time gave noisy results, and the longer times tended to 
mask the degree of agreement between the two field orientations. 

Examples of the results are shown in Figure 2. Solar wind speeds 
ranged from 340 to 425 km/s during these intervals. The heavy traces give 
the time variations of the angle between the total field vectors, and the 
thin traces give the variations between the projections of the vectors in 
the y-z plane. The lag times of the minimum angles between the total and 
projected fields at the two spacecraft are labeled and T^ ^ , 
respectively. 

The five panels in Figure 2 illustrate five types of variations. In 
the top two panels the interval of lag time over which the angles are 
minimal is relatively short and well-defined, especially in the top panel. 
Plots of this type meet the commonly held expectation that magnetic 
features observed far upstream are convected downstream with little 
change, and that their scale size in the outflowing direction is small 
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relative to the interspacecraft distance. The third panel shows evidence 
of large-scale wove structure with a period of 30 min and a wavelength 
of - 100 Rg. In cases like this the angle between the fields at the two 
spacecraft reaches minimal values at several different lag times. In the 
fourth panel the IMF has essentially the same orientation at the two 
spacecraft throughout the two-hour period. Thus the scale size of the 
magnetic feature being observed must be considerably larger than the 
Interspacecraft distance. The fifth panel shows what appears to be a tran- 
sition from a period like that in the third panel to a time when the fields 
at the two spacecraft are completely uncorrelated. The lack of correlation 
may be owing either to temporal changes or to scale sizes perpendicular to 
the flow direction which are less than the interspacecraft distance in the 
y-z plane. In the total set of cases, variations of the types shown in 
panels two, three, and four are about equally common. Only 37 . of the cases 
show the poor correspondence of the fifth panel, and only 6% show the excel- 
lent definition of the first pattern. 

The vertical lines in Figure 2 mark two predicted lag times. The solid 
line marks the simple convection lag time T conv = Ax/V Ce.g., King, 19833, 
where Ax is the difference in the x coordinates of the two spacecraft posi- 
lons, and V is the solar wind speed. In using T conv as a predicted lag 
time, one assumes that surfaces of constant phase in a magnetic feature lie 
perpendicular to the radial solar wind flow direction. The dashed line 
marks the lag time based on the IMF orientation, Tijjp. It is the same as the 
corotation lag time Ce.g., King, 19833 except that it is calculated from the 
measured IMF orientation at ISEE 3 rather than the assumed Archimedean spi- 
ral orientation. 

Figure 3 illustrates how Ti^p is determined. The plane of the figure 


contains the ISEE 3 position and is parallel to the x-y plane of the geo- 
centric interplanetary medium <GIPM) coordinate system CBieber and Stone, 
19793. Spacecraft coordinates and magnetic field vectors are transformed 
from the GSE to GIPM coordinate system by a rotation about their common x 
axis. The angle of rotation is determined by the requirements that the IMF 
lie parallel to the x-ygiPM Plane <B z »gxPH 5 0) and that its x and y compo- 
nents have opposite signs. The sign requirement causes the IHF to be 
directed toward or away from the earth from the "dawn" <+x,-y) quadrant. 
An IMF vector f is drawn through the ISEE 3 position in Figure 3. If it is 
carried to a projected ISEE 1 position at speed V, then the time it takes to 
reach that position depends upon the y coordinates of the spacecraft posi- 
tions as well as the x coordinates. Three different projected ISEE 1 posi- 
tions are shown, at x distances Lx, 1*2, and L3 from 5. Since Tijjp » L/V, 
and Lx < L2 < L3, it follows that TxHpv.x < Txhf -2 < Tjjjp-3. In the case 
of position 2, the two spacecraft have the same y coordinate, so that 
L2 = Ax, and Tjnp_2 = T conv * In general the ISEE 1 spacecraft lies either 
above or below the plane of the figure. Thus in using Txmf 83 a predicted 
lag time, one assumes that surfaces of constant phase in a magnetic feature 
lie perpendicular to the x-yQipn plane as well as parallel to 3. 

In the top panel of Figure 2, Tihf is nearly identical to T conv . Near 
coincidence of Txmf an< l T conv can occur either when the yGIPM distance 
between the spacecraft is small, as in the center position in Figure 3, or 
more commonly when the x component of the IMF is small. What is remarkable 
about the predicted lag times in the top panel is that they also coincide 
with the lag times of minimum angle, T m x n . In the second panel the pre- 
dicted lag times are widely separated: Tjjjp differs from T conv by 31 min. 
The former clearly is the better predictor in this case, since it lies 
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wlthln the veil of minimum angle between the two spaccraft IMF orientations 
In the remaining panels, Tj^p n T conv , ns in the top panel, but they do not 
coincide with the T m i n values. The predicted logs in the third panel indi- 
cate which of tb> several minimum angle wellB in the wave structure corre- 
sponds to the same solar wind parcel passing the two spacecraft. Inter- 
estingly, it is not the well containing either of the T m i n values. Although 
T C onv ar, d Tihf are within about 5 min of what appears to be the true lag 
time at the nearest local minimum, the differences between T conv or and 
T K ln are much larger and contribute unwarranted spread to the statistical 
distributions of lag time differences shown later. In the fourth and fifth 
panels the T n i n values are relatively meaningless. Even though the IMF 
orientations at the two spacecraft ore well-correlated in the fourth panel, 
the lack of time variation of the angle between them affords no Information 
about lag times. Cases like this also contribute large spread to the time 
lag difference histograms but sharpen the histograms of angular differences 
at predicted lag times. The fifth panel has large time variations, but lag 
time comparisons are not useful here either because of the Increasingly poor 
correlation between the two IMF orientations throughout the interval. 

Statistical Distributions of angular Differences . Figure 4 shows histo- 
grams of the angle between the IMF at ISEE 1 and ISEE 3 at no lag time T„ 
<top panels), at the convection lag time T conv (middle panels), and at the 
IMF orientation lag time Tj^p (bottom panels). The left panels give the 
angles between the total IMF vectors, and the right panels give the angles 
between their y-z projections. The angles below which SOZ (medians) and 80JJ 
of the cases fall are listed in Table 1 along with the mean angle for each 
histogram. 
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The gentle peaks on the distributions of angles at T 0 indicate that the 
fields remain somewhat correlated even when no transport time is token into 
account. Similar results were found by King C 1983]. Half of the cases 
have angles less than ~ 35°, and 8031 less than ~ 6S°, both for the total 
and the projected IMF vectors. However, these statistics are improved 
considerably at the predicted lag times. At T conv , half of the coses have 
angles less than 20°, and 802 less than 38°. At Tiny, the angles are 

reduced only by ~ 1°. 

Additional histograms are shown In Figure S. At the top are the 

angular distributions at the lag time T m i n . They give the smallest 
angular difference between the IMF at the two spacecraft within each 

two-hour data segment. Thus they represent an upper limit on how well 
the IMF near the earth can be predicted by using time-lagged measure- 
ments from a librati on-point-orbiting spacecraft. Half of the cases 
have angles less than - 10°, and 80Z less than - 15° (see Table 1). A 
comparison between these histograms for T m i n and the histograms for 

Tconv ln ^e middle panel, r produced from Figure 4, indicates the 
degree of uncertainty introduced by using a predicted lag time. The 

histograms in the bottom panel are a quantitative measure of that 

uncertainty. They give the distribution of the difference between the 

minimum angle and the angle at the predicted lag time. (Note that 

these difference distributions are not the result of a simple subtraction 
of the above distributions.) Table 1 indicates that the median angles at 
T CO nv are ~ 8° larger than the median minimum angles. 

In order to test the sensitivity of the angle at T conv to lag time, a 
histogram of the smallest angle within S min of T conv was constructed (not 
shown). The mean of this distribution is 37°, 6° less than the mean of the 






distribution for T conv . This value translates roughly into an average slope 
of 35° per half hour of lag time for the sides of the veils of minimum 
angle in the lag time versus angle plots, as in Figure 2. 

Figure 6 shows normalized angular distributions of cases selected 
according to whether the interspacecraft distance D in the y-z plane was 
less than SS % or greater than 90 %. The top, middle, and bottom panels 
give the angles between the total INF vectors at the two spacecraft at lag 
times T conV f Tjjif, and T m i n , respectively. It is clear that the angles ore 
smaller for smaller D, in agreement with Crooker et al. [19821. The lower 
half of Table 1 lists the medians, means, and 80Z levels of the distribu- 
tions for the two categories. Overall, the values for small D are ~ 30Z 
better than those for large D, with the most improvement <« 3SZ) for the 
T conv distributions. Table 1 also lists the medians, means, and 80 Z levels 
for distributions of the angles in the y-z plane (not shown). The improve- 
ment in these values for smaller D is less, ~ 1S-20Z. 

Time Lag Statistics . The average time for a magnetic feature to travel 
from I SEE 3 to ISEE 1 was ~ SS min during the period studied. This value 
was determined from histograms (not shown) of the lag time T m i n , when the 
angle between the IMF vectors at the two spacecraft reached its minimum 
value. The median and modal values of these distributions for both the 
full and projected vectors are S3 min. 

Histograms of the difference beween T m ^ n and the predicted lag times 
T C onv and Tijjp are shown in Figure 7. They are approximately symmetric 
about zero and peak there. Thus the predicted lags are neither systemetic- 
ally too long nor too short. The peaks for T m i n - Tj^p rise more sharply 
than those for T m j n - T conv , indicating that Tjjjp is a somewhat better pre- 
diction. For example, in the total vector histograms on the left, 48 Z of 
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the T m in - Tijjp values lie within ± 10 Bin of zero, compared to 402 of the 
^inin ” T conv values. 

These results suggest that IMF features ten# to conform to the geometry 
in Figure 3, for which Tjjjp was designed. In order to test the strength of 
this tendency, the T m i n - T conv histograms in the top half of Figure 7 were 
separated according to whether L is less than or greater than Ax. If L < 
Ax, as is Li in Figure 3, then the observed lag time T m ^ n should be less 
than T conv , since the solar wind has a shorter distance to travel to bring 
the magnetic feature S to ISEE 1 than to travel the full x distance between 
the spacecraft. If L > Ax, as is Lg In Figure 3, then T n i n should be 
greater than T conv . 

The results of this test are shown in Figure 8. In order to reduce the 
noise level, cases with IL-Axl < 10Rg were eliminated. The top panel is 
for the total field vector differences, and the bottom for the y-z projec- 
tions. It is clear in both panels that the distribution of cases with L < 
Ax peaks on the negative side of the T m ^ n - T conv axis, and vice versa, in 
agreement with the prediction. 


Discussion and Conclusions 

The main conclusion to be drawn from this study is that the IMF orien- 
tation near the sunward libration point, lagged simply by the time 
required for the solar wind to convect to the earth, is a reasonably good 
predictor of IMF orientation near the earth. Specifically, there is a 502 
(802> probability that a ten-minute averaged predicted IMF vector lies 
within 20° (38°) of its true direction. . The probabilities are essentially 
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the same for the y-z vector projections. 

A second conclusion Is that a lag time T^, based on the assumption that 
magnetic features tend to be aligned with the observed field direction, Is 
clearly more accurate than the simple convection lag time T CO nv- However, 
on a statistical basis, the difference in the two predicted lag times is 
not large enough to improve significantly the distribution of angular 
differences between the IMF vectors at the two spacecraft. ThU 3 , since 
T conv ls nuch nore convenient to calculate than Tngr, the IMF orientation 
at T conv remains the predictor of choice. 

The uncertointy In the predicted IMF orientation has two components. 
The first results from uncertainty in the log times The angle between the 
llbratlon point IMF vector and the near-ea**th vector ls always larger at 

Tconv thQn at T mln , unless T conv = T B i n , since by definition T Bln ls the 

lag time when the angle reaches its minimum. There is a SOZ <0OZ> 

probability that this angular difference for a ten-minute averaged vector 
is less than 0° (20°i. The second component of uncertainty results from 
temporal changes in IMF orientation as the solar wind convects from the 
llbratlon point to the earth, and also from spatial scales of magnetic 
features which are smaller than the y-z projection of the distance from 
the earth to the libration-point-orbiting spacecraft. It is these 

variations which prevent the angle between IMF orientations at the two 
spacecraft from reaching zero at T mln . There is a SOZ (80Z) probability 
that a ten-minute averaged IMF vector at the llbratlon point will lie 
within 11° (19°) of the orientation of a ten-minute averaged vector near 
the earth during the following two-hour period. Thus the two components 
of uncertainty are about equal and account for the total uncertainty In 
the predicted IMF orientation given by the probabilities above. 


The final conclusion concerns the Improvement In the predicted IMF 
orientation vhlch may be obtained by restricting the data to those times 
when the libratlon point orbit Is within some distance D of the earth-sun 
line. As mentioned in the introduction, Crooker et al. C 19023 found that 
the correlation between fields at ISEE 3 and ISEE 1 falls sharply for D > 
90 Rj? . The results here indicate that the angles between the fields at the 
two spacecraft are about 30% smaller for D < 55 R£ than for D > 90 R E . 
However, the statistical improvement gained by simply restricting cases to 
D < 90 Re is considerably less, about 10%. Specifically, there is a 50% 
<80%) probability that the predicted IMF vector is within 18° (33° > of its 
true direction In the restricted data set, compared to within 20° (38°) for 
the total data set. 

Although the subject of scale sizes of magnetic features Is covered 
elsewhere Ce.g., Crooker et al., 19821, additional information on scale size 
along the earth-sun line Is obtained from this study. The histograms of 
angles between the IMF at the two spacecraft for no time lag (top of Figure 
4) indicate thot features often span the x distance between the spacecraft, 
~ 200 R e . Further, scale sizes greater than 400 R E along the earth-sun line 
(x-axis) must have been present about a third of the time during the period 
studied. Roughly a third of the cases had essentially the same IMF orienta- 
tion at the two spacecraft over the length of the two-hour data segments 
analyzed, as In the fourth panel of Figure 2, and the solar wind normally 
covers a distance of more than 400 Ft E during two hours. 

The period of time covered by the study was near the maximum phase of 
the solar cycle, when large scale coronal mass ejections occurred at a rate 
of - 1/day CHundhausen et al., 19841. These ejections may account for the 
large scale sizes deduced from our study. Since the scale size of their 
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magnetic features Is Inferred to be ~ S000 Rg Ce.g., Klein and Burlnga, 
1982; Crooker, 19831, only about 1/3 of the ejections would be required to 
pass In the earth’s vicinity to account for the frequency of the signatures 
in the fourth panel of Figure 2. 
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TABLE 1. Characteristics of Distributions of 
Angular Differences Between IMF Vectors at ISEE 1 
and ISEE 3 


Total Vector / y-z Projection 

Las Time Mean Median 80% Level 

ALL CASES (885) 


T 0 

36°/34° 

^conv 

20°/18° 

t imf 

19°/16° 

T 

min 

ll°/0° 

T - T 

conv min 

8°/9° 

X I1IF " T min 

7°/8° 


CASES 

^conv 

17°/1S° 

t imf 

16°/ 14° 

T 

min 

9°/7° 


CASES 

^conv 

26°/20° 

t imf 

23°/17° 

T 

min 

130/80 


43°/46° 

63°/69° 

26°/27° 

38°/38° 

2S°/2S° 

37°/37° 

13°/10° 

19°/14° 

13°/1S° 

29°/2S° 

12°/1S° 

18°/22° 

WITH D < SS 

Rj. (414) 

22°/24° 

31°/34° 

22°/22° 

30°/32° 

11°/10° 

17°/14° 

WITH D > 90 

(221) 

32°/31° 

49°/48° 

29°/28° 

45°/40° 

1S°/11° 

21°/16° 
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Fig. 1. ISEE 3 orbit In GSE coordinates from July 25 to Nov. 27, 1979. 
The shaded region in the upper panel Indicates the simultaneous coverage 
of ISEE 1. 


Fi g . 2. Examples of variations In the angle between the IMF at ISEE 1 and 
ISEE 3 as a function of lag time. The solid traces give the angle between 
the total vectors, and the dash-dot traces between the projections of the 
fields in the y-z plane. The vertical solid and dashed lines mark the 

predicted lag times T and T.,_, respectively. Arrow labeled T . mark 

the lag times of minimum angle between the total vectors (B) and between 
the y-z projections <yz). 


Fig. 3. Schematic diagram illustrating the geometry for calculating Tj^. 
The plane of the diagram is parallel to the xTqjpjj plane and contains the 
ISEE 3 position, through which an IMF vector 3 is draw. Three possible 
projected ISEE 1 positions 1, 2, and 3 are shown (not to scale). The vector 
& convects to positions 1, 2, and 3 over distances 1^, L 2 , and L g at the 
solar wind speed V. Since < L 2 < L gf the transit time Tj^ is shortest 
for position 1 and longest for position 3. At position 2, = Ax/V = 
T C onv' where Ax is the difference in x coordinates of the ISEE 1 and ISEE 3 
positions. 


Fig. 4. Histograms of the angle between the IMF at ISEE 1 and ISEE 3 at 

zero lag time T and at the predicted lag times T and T TUrt . Angles 
o conv IMF 

between the total vectors B are on the left, and angles between their y-z 
projections are on the right. 


Fig. S. Histograms of the angle between the IMF B vectors and y-z 

projections at ISEE 1 and ISEE 3 at T , from Figure 4, and at T , the 

conv win 

lag time when the angle reaches its minimum in each two-hour data segment. 

In the bottom row are histograms of the difference between the angle at 

T and the minimum angle, at T , 

conv min. 


Fig. G. Comparison of histograms of the angle between the IMF B vectors at 

ISEE 1 and ISEE 3 at T . T tm1? , and T . for D < 55 R-, (414 caseB) and D 

conv IMF min b 

> 90 Rj, (221 cases), where D is the interspacecraft distance in the y-z 
plane. The histograms are normalized to percentages to facilitate 
comparison. 


Fig. 7. Histograms of the differences between the lag time at minimum 

angle T . and the predicted lag times T and T t1 _, for the IMF B vectors 
min conv xnr 

and y-z projections. 

Fig. 8. Comparison of histograms of T Bln - ^ conv for < *** cases) 

and L > Ax (194 cases). These conditions are illustrated at positions 1 and 

3, respectively, in Figure 3. Cases with IL - Ax I < 10 R^ were eliminated 

to reduce noise. As predicted, cases with L < Ax have shorter lag times 

T . than the predicted lag times T , and vice versa. 

min conv 
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Figure 6 













